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Abstract 

We explicitly compute the tree-level on-shell four-graviton amplitudes in four, five and six dimensions for 
local and weakly nonlocal gravitational theories that are quadratic in both, the Ricci and scalar curvature 
with form factors of the d’Alembertian operator inserted between. More specifically we are interested in 
renormalizable, super-renormalizable or finite theories. The scattering amplitudes for these theories turn 
out to be the same as the ones of Einstein gravity regardless of the explicit form of the form factors. As a 
special case the four-graviton scattering amplitudes in Weyl conformal gravity are identically zero. Using 
a held redefinition, we prove that the outcome is correct for any number of external gravitons (on-shell 
n —point functions) and in any dimension for a large class of theories. However, when an operator quadratic 
in the Ricmann tensor is added in any dimension (with the exception of the Gauss-Bonnet term in four 
dimensions) the result is completely altered, and the scattering amplitudes depend on all the form factors 
introduced in the action. 
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1. Introduction 

Scattering amplitudes, in particular gauge invariant on-shell ones, are the observables of major interest in 
particle physics experiments. The possibility of getting very precise and well testable predictions for them is 
surely one of the biggest achievements of perturbative quantum field theory. From a more theoretical point 
of view, they play a fundamental role in constraining the higher derivative terms that are expected to arise 
at quantum level as a consequence of the presence of divergences in loop diagrams in non-renormalizable 
theories [l|]. It is well known that, even for tree-level diagrams, the intermediate steps of computation show 
a remarkable tangle that often disappears in the final physical amplitudes @. This is a consequence of the 
ambiguity (and also the richness) intrinsic to the Lagrangian off-shell formalism, which is invariant under 
gauge transformations and field redefinitions. Recently there has been a wide interest in developing powerful 
methods for on-shell scattering amplitudes, which have led to both improving calculation techniques and 
gaining new insights into the underlying mathematical structures (see [4[ and references therein). 

On the other hand, since it is quite difficult to propose good and unambiguous quantities which could 
play the role of observables in pure quantum gravity, we concentrate our attention on hypothetical exper¬ 
iments of graviton scattering described in the perturbative framework of quantum field theory around flat 
Minkowski spacetime. Due to the kinematical constraints the first non-vanishing amplitude (for on-shell 
massless gravitons) must describe a scattering process involving four particles. One example is the scat¬ 
tering of two gravitons into two, which is the topic of this paper. Here due to tremendous complications 
present at loop levels in quantum gravity, we consider only tree-level amplitudes (only in the Born approxi¬ 
mation). Physically these amplitudes correspond to the classical processes of monochromatic gravitational 
wave scattering happening in the empty space. Of course such effects are predicted by Einstein general 
relativity, which is a nonlinear theory of gravitational interactions, but non-linearities are very little effects 
due to the smallness of the ratio between the typical graviton energy and the Planck mass. Nonetheless, the 
amplitudes are indeed very useful probes of quantum gravitational physics and are the first must-be-taken 
step to quantum gravity phenomenology. 

The first attempt to compute tree-level scattering amplitudes was undertaken in the simplest quantum 
gravity theory, namely in the quantum version of Einstein gravity. There it was found (3, (§.] that the 
final results are quite simple and can be derived even without performing detailed computations using 
Feynman diagram methods [BJ. Due to dimensional reasons (Gjv has energy dimension —2 in four dimensions) 
the amplitudes showed behavior that grow unboundedly with increasing energy of the scattering process. 
Precisely they grow as E 2 . This observation led to the conclusion that naive tree-level unitarity bound on 
the scattering amplitudes is violated. Further, it was derived that either Einstein-Hilbert quantum gravity 
becomes non-perturbative around Planck scale or the theory needs an ultraviolet (UV) completion and 
can be viewed merely as a low energy effective field theory of quantum gravity. We keep this message in 
mind and we study the scattering processes in a class of theories which can be without problems in the 
UV regime, but at the same time these theories remain always in the perturbative regime. Indeed, higher 
derivative super-renormalizable theories are asymptotically free and the unitarity bound is not violated. 
This is transparent in the prototype asymptotically free theory proposed by Stelle in 1977 @-[s]. The energy 
scaling of the scattering amplitudes is proportional to E 4 or E 2 , but at a closer inspection it can be correct 
only up to the Planck scale because at very high energy the interaction between gravitons becomes very 
weak and they travel almost as free particles as a mere consequence of asymptotic freedom. Therefore, 
the scattering amplitudes presented in this paper are meaningful up to the Planck scale (or the new scale 
introduced in the theory), while a careful analysis is required around the Planck mass. 

Moreover, despite the remarkable result of our computation, which states that amplitudes in four dimen¬ 
sions are the same as in Einstein gravity, we could be able to determine other phenomenological consequences 
of the graviton scattering processes. One of such interesting implications concerns the production of micro 
black holes as resonances in the elastic scattering of gravitons |9|, 110|. Typically, as it is done in high energy 
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physics, from a broadening of the cross section we would be able to read out the half-width and therefore 
determine the life-time of such objects. This would shed a new light onto the interpretation of Hawking 
evaporation process in the fully-fledged (opposed to semi-classical approach) quantum field theory of gravi¬ 
tational interactions. Other implications for particle physics and its unification with gravity would require 
the coupling of matter to gravity, which we do not attempt to make here. 

As the last piece of our original motivation we want to mention that the results of the amplitude 
computation might have served well for the determination of form factors, which are used in nonlocal 
theories (for the special case of constant form factors we would be able to determine values of the couplings 
in front of terms quadratic in curvatures in the Stelle gravity). But since in four spacetime dimensions 
we find no dependence on interpolating functions (form factors) with the exception of terms quadratic in 
Riemann tensors, this cannot be done based only on tree-level results for the minimal super-renormalizable 
or finite theory. The situation is different in higher dimensions as we point out explicitly in the second 
part of the paper and this is one of the new results of this paper. In the higher dimensional setup, thanks 
to these amplitudes, we are able to determine some characteristic features of the form factors. For fully 
unambiguous determination of the form factors in the minimal nonlocal super-renormalizable theory we will 
need the full one-loop result, which is right now beyond our capabilities. It can be quite pleasing that in 
four physical dimensions tree-level amplitudes do not depend on the form factors, because then we can view 
the presence of the form factors as the clever way of parametrizing our ignorance. However, as we see on 
the tree-level nothing depends on it and actually we are not forced to determine the form factor there. This 
changes one of the biggest drawback of the theory (ambiguity related to the choice of the form factor) into 
a virtue (now tree-level scattering amplitudes are independent of it). This further may lead to seeing the 
form factors as something that does not have very important physical meaning. The form factor may be 
used only as a way to interpolate between infrared (IR) and UV behavior of the theory in a covariant way. 
Such ideas are similar to the ideas of the smooth covariant cutoff developed in the field of functional RG. 
We believe that the true significance of the form factors used in nonlocal theories will be revealed only at 
the loop level. This is confirmed by partial computations of the beta functions of running couplings, which 
indeed do depend on the asymptotic form of the form factors used. 

All of this has drawn our attention to investigations, with somewhat more traditional techniques, of 
what kind of physical information can be obtained from tree-level on-shell amplitudes for a class of weakly 


nonlocal theories of gravity, which have been the subject of recent studies 11 m . These theories turn out 
to be a particularly convenient theoretical framework to perform quantum gravity computations, because 
they have been proven to be ghost-free and super-renormalizable or finite at quantum level. We work in 
the quantum field theory framework and we assume as our guiding principle the “validity of perturbative 
expansion” [21[. Moreover, the following postulates are required: (i) spacetime diffeomorphism covariance; 


(ii) weak nonlocality (or quasi-polynomiality); (iii) unitarity (the spectrum is tachyon- and ghost-free); 
super-renormalizability or finiteness. The main difference with perturbative quantum Einstein gravity lies in 
the second requirement, which makes possible to achieve unitarity and renormalizability at the same time. 
In this way a fully consistent quantum gravitational theory has been constructed free of any divergences. 
However, the theory is not unique and all our ignorance is encoded in a form factor (entire function) with 
very specific asymptotic limits in the ultraviolet and in the infrared fixed in such a way to have a convergent 
quantum field theory. At classical level there are evidences that we are dealing with a “singularity-free” 
gravitational theory in the case of physical matter [22l - f29l |. However, we know that Einstein spaces are exact 
vacuum solutions of the weakly nonlocal theory, including the singular Schwarzschild spacetime 0,E1. 

The main findings of this paper consist of new results for tree-level scattering amplitudes in Stelle gravity 
and its analytic nonlocal extensions. We find that the situation in quadratic gravity in four dimensions is 
exactly the same as in Einstein gravity, while for the Weyl conformal theory the four-graviton amplitudes 
are identically zero. In dimension higher than four and when the terms quadratic in Riemann tensors are 
included the outcome in Stelle gravity differs from the standard Einstein theory. If in four dimensions (or 
extra dimensions) we allow for analytic nonlocal extensions involving the Riemann tensor, then again we find 
a different outcome. We filially explain our results by proving a general theorem about on-shell n —graviton 
scattering amplitudes and by explicitly considering diagrams from terms that cannot be redefined. We 
point out the possibility of determining the form factors, when these theories are viewed as fundamental, 
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and comparing them with the experimental results about graviton scattering amplitudes. This may be taken 
as a first step to verify this class of theories. 

In this paper we mainly concentrate on the calculation of tree-level four-graviton scattering amplitudes. 
In the remaining part of this introduction we remind the reader the construction of the weakly nonlocal 
quantum gravity. Later, in the second section, we start with some kinematical considerations and we find the 
propagator and compact expressions for vertices. We first perform the amplitude computation in quadratic 
Stelle theory [§] in section three and then in the most general weakly nonlocal super-renormalizable theory 
quadratic in the Ricci and scalar curvature tensors in section four. Moreover, in the fifth section we interpret 
the technical results as a consequence of a general theorem based on Anselmi’s field redefinition theorem 
: 32, 331], and [MS- Finally, we discuss the generalizations of our results (also some expectations beyond 
tree-level) and we give conclusions. We supplement the paper by two appendices, where we put more 
technical details about the propagator and variations of curvature invariants. 


1.1. Weakly nonlocal gravity 

The general U-dimensional theory weakly nonlocal and quadratic in Riemann, Ricci and scalar curvature 
reads [0, 0,0, ES E3 0, SBl , 


C g = —2 k d 2 \/—g [R + R- 7 o (□)!?. + Ric 72 (D)Ric + Riem 74 (D)Riem + V] . 


(1) 


The theories above consist of a weakly nonlocal kinetic operator and a local curvature potential V. Therefore 
they are quite general, more general thing would be only to allow the potential to be nonlocal too. The 
local potential V for a gravitational theory, which is however sufficient, is made up of the following three 
sets of operators 


N+2 j 


7+N+l j 


E E E 43 (y 2U ~ k) n k ) i + E E E 43 (v 2(j " 


k) n i 


j =3 k =3 


j— N+3 k —3 


7+N+2 

fE E*M (v 2( ^+ N+2 - fc) ^ fc ) i , ( 2 ) 


fc=3 


where operators in the third set are called killers, because they are crucial in making the theory finite in 
any dimension. They are used to kill the beta functions. In @ Lorentz indices and tensorial structure 
have been neglectcclJ. For more details about notation in the potential V we refer the interested reader to 
original papers fl2l. Tl3i] . Moreover in O □ = v is the covariant box operator, and the functions 

7 f(D) (£ = 0,2,4) (form factors) will be shortly defined. When 74 (D) = 0 minimal unitarity requires 


7o(D) 


72 (D) - 1 

2 ~ 2D 


( 3 ) 


where H(z) is an entire function on the complex plane. Given the conditions Q, in the perturbative 
spectrum of the theory we only have the massless transverse graviton, whose propagator is modified only in 
a multiplicative way from the Einstein theory. In 0 A is a fundamental invariant mass scale in our theory 
and we define z = —Da = —D/A 2 . 

A universal exponential form factor exp H (z) compatible with the guidin g pr inciples of quantum field 
theory and with the requirement of proper IR limit of the theory H( 0) = 0 is [171 ]: 


V~ 1 (z) = e = e t[ r ( 0 ’PG) 2 )+'i' E + 1 °s(p( z ) 2 )] 


= e a ~^~ 


a e 

2 p{z) 2 



( 4 ) 


where the last equality is correct only on the real axis. 7 e ~ 0.577216 is the Euler-Mascheroni constant and 
r(0, z) = /„ +0 ° dte^/t is the incomplete gamma function with its first argument vanishing. Now we define 
capital N by the following function of the spacetime dimension D: 2N + 4 = D in even and 2N + 4 = D + 1 


1 Definitions — The metric tensor f] jlu has signature (-(-■••-(-) and the curvature tensors are defined as follows: R^ pc7 = 

—dcrVvp + ..., Rju/ = R p pv, R = R P v- With symbol 7Z we generally denote one of the above curvature tensors. 


4 

















in odd dimensions. Besides integer N in our theory we have also another integer number 7 , which measures 
how far the theory is from the minimal renormalizable one in given spacetime dimension. The polynomial 
p(z) in (01) is of degree 7 + N + 1 and such that p( 0) = 0, which gives the correct low energy limit of our 
theory (coinciding with Einstein gravity). The entire function has asymptotic UV behavior given by the 
polynomial p(z) a in conical regions around the real axis. However in most applications we will choose a 
monomial behavior ~“(t+n+i)^ w } iere a j s another positive integer number. Then these conical regions are 
with an opening angle 0 = 7 r/( 4(7 + N + 1)). For 7 = N = 0 we have the maximal angle 0 = 7 r /4 for all a. 


1.2. Propagator, tree-level unitarity and power counting 

We shortly discuss the propagator in a theory ©• Splitting the spacetime metric into the Minkowski 
background plus fluctuations and after performing a gauge fixing, we can invert the quadratic in fluctuation 
fields kinetic operator to finally get the two-point function in Fourier space j50| (see section (12.21) and 
|Appendix A| for more details), 

V{k 2 / A 2 ) P oe 


o~ = - 


k 2 


P 2 - 


D -2 


(5) 


The indices for the operator 0 _1 and the projectors [50 : , 53] {P 09 , P 2 } as well as gauge-dependent terms in 
(|5|) have been omitted. The tensorial structure in © is the same of Einstein gravity, but the multiplicative 
form factor V{k 2 / A 2 ) in Fourier space makes the theory strongly ultraviolet (UV) convergent without the 
need to modify the spectrum or introducing ghost instabilities. The theory propagates only massless spin 
two particle (graviton with two helicities), however the UV behavior of the propagator is modified. This 
fact about the tree-level spectrum ensures us that we are dealing with a perturbativelv unitary theory. 

11 nay, 0." 


We now review the power counting analysis of the quantum divergences [11| 


In the high 

energy regime, and in even dimension the above propagator m in momentum space schematically scales as 


0~\k) 


1 


k^J+ D 


in the UV. 


( 6 ) 


The interaction vertices can be collected in different sets, that may or may not involve the entire functions 
exp H(z). However, to find a bound on the quantum divergences it is sufficient to concentrate on the leading 
operators in the UV regime. These operators scale as the propagator giving the following upper bounds on 
the superficial degree of divergence of any graph H, 46], 


w(G) = DL + (V - /)(27 + D). 

We rewrite above in a more convenient form as 
oj{G) = D — 2 7 (L — 1). 


(7) 


( 8 ) 


In ([5]), we used the topological relation between the numbers of vertices V, internal lines I and the number 
of loops L: I = V + L — 1. Thus, if 7 > D[2_. in the theory only 1-loop divergences survive. Therefore, the 


theory is super-renormalizable 0, 0, ]23, 45:] and only a finite number of operators of mass dimension up 
to M D has to be included in the action in even dimension. For odd dimension, if 7 > (D — l)/2, then the 
theory is completely without divergences and hence automatically finite. 


1.3. Super-renormalizable & finite gravitational theories in D = 4 

The main reason to introduce a potential V in the action m is to make the theory finite at quantum 
level. It is easy to see that it is always possible to choose the non-running coefficients Sk,i in m to make all 
the beta functions vanish. We consider the simplest case of a monomial asymptotic behavior for the form 
factor, namely: p 7 + n+i(z) = £ 7+N+1 . For this particular choice of the form factor in the large 2 : limit, the 
analysis of the previous subsection shows that only divergences, which are to be renormalized by terms with 
D derivatives (like 'RP^ 2 /e in dimensional regularization (DIMREG)), are generated at one-loop. However, 
the killer operators in the last set of operators in m give contributions to the beta functions of the theory 
linear in the front parameters Sk,i■ The latter ones can be fixed in such a way to make the theory finite 
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in any even dimension. The tensorial structure of killers must reflect the structure of terms, which are 
renormalized in the original theory (terms with D derivatives). 

In D = 4 (N = 0) the whole situation is simple to describe, because we only need two killer operators. 
The highest derivative terms in the kinetic part of the action come from the form factor and are of the type 
1Z UPTZ [48l450l ] . The minimal choice for a finite and unitary theory of quantum gravity in four dimensions 
may consist of terms with 7 = 3 (and a = 1) in the kinetic part. This alone leads to one-loop super- 
renormalizable quantum nonlocal gravity. For simplicity we introduce only two quartic killers and no cubic 
in curvature operators. This is sufficient to make vanish all two beta functions for R 2 and R 2 IV operators. 
The simplest Lagrangian with the Tomboulis type of form factor (j4j) may be the following, 




R + R 


■fAlS 


□ 


R! lv - R 


s if(—n A ) _ \ 
2D 


R + Sl R 2 a(R 2 ) + s 2 R^R^a(R P aR p,T ) , (9) 


where p{z) = z 4 , Si = —27r 2 W2(ci + C2)/3, S2 = 87t 2 W2C2 and w2 = A -8 exp(7s/2) [l2j]. Here c\ and C2 
are two constants independent on u> 2 , that have to be determined from the calculation of the contributions 
to the beta functions from terms quadratic in curvature and dominant at high energies. A more general 
Lagrangian can have a finite number of other local terms, but still finiteness of the theory can be obtained 
exactly in the same way. 

We recently proposed, following [20.], another class of weakly nonlocal possibly finite theories, which 
are constructed entirely from kinetic terms (only weakly nonlocal operators quadratic in curvature appear), 
without local or nonlocal gravitational potential V cubic in curvature or higher. The simplest nonlocal four¬ 
dimensional theory we can write to achieve finiteness is ex actly the one given in ©, but without assuming 
the relation © (unitarity is anyhow achieved as proven in [ 14, 0 ), 

( 10 ) 


£ = -2k 


7 ’<(□) = 


fV=g[R + R^o(a)R + r^ 2 {u)r pv + R^ lA {n)R^) . 

e H t (- Da) _ 1 


□ 


We assume that all three form factors (for £ = 0,2 and 4) have the Tomboulis form (0]) with the same degree 
of UV polynomial 7 + I. The operators in m can be written equivalently in other bases using Weyl tensors 
C fJjU pa r )wi^)C ll '' pCT or generalized Gauss-Bonnet Lagrangian. The finite theory would boil down to some 
relations between UV polynomials for each form factor. At the moment we cannot definitely assert that the 
theory is finite, because the beta functions for the operators R 2 and R 2 ^ v are quadratic in three parameters 
defining the asymptotic UV behaviors of the entire functions Hg. This theory and its generalizations are 
currently under careful investigation and the results will be published in a future paper. Nevertheless, 
looking at the beta functions reported by Kuz’min j20], it turns out that for 7 > 4 we can always find a 
solution for 


P R 2 = 0 or P R 2 ^ = 0 , 

hence in this class of theories only one beta function could be non-zero. 


( 11 ) 


2. Four-graviton scattering amplitudes in higher derivative gravity 

We start by studying the higher derivative gravity theory proposed and extensively studied by Stelle in 
1977. The action is quadratic in curvature and it is the limiting case of m for the particular case of constant 
form factors and zero potential, namely: 

7o(D) = const = 70 , 72(D) = const = 72 , 74(D) = const = 74 , V = 0 . ( 12 ) 

Therefore, the higher derivative Stclle’s gravitational action reads 

£ g = —2 k d 2 Jd D x \/~^g [i? + 70 R 2 + 72 R 2 ^ + liR^upa] ■ (13) 

In this section, before carrying out the computation of the on-shell four-graviton scattering amplitudes, we 
derive the propagator and vertices and discuss some general properties of helicity amplitudes. 
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2.1. Helicity amplitudes 

For the external on-shell gravitons we assume the physical conditions in transverse traceless gauge 

^ V = K = 0 ’ ( 14 ) 

which will turn out to be very convenient in order to simplify the algebra all along the computation. In this 
gauge i the polarization tensors for gravitons with helicities ±2 in four dimensions are related to those of 
photons with helicities ±1 by tensor product 

e^(p,±2) = e^{p,±l)e„{p,±l), (15) 

where the polarization vectors satisfy 

<+(p, A)p M = °> <+0, A)e M (p, A) = 0 , 

<V(p,-A) = e*(p,A), e M (p, A)e M (p, —A) = 1. (16) 

Furthermore, we require the polarization vectors to form a complete basis for a representation of the SO( 2) 
group of rotations in the transverse directions which leave p^ invariant. In the axial gauge we introduce the 
auxiliary vector q^ ])[ p ^ such that q ■ p ^ 0, and the sum on the polarizations reads 


^2 A )e*(p, A) = v - 

A=±l 


P^v + T+9 M 
p-q 


(17) 


In particular, taking momenta in spherical coordinates (0, cp) to be 
TV = {po > Po sin @ cos <Pi Po sin 6 sin <p, po cos 9) , 


(18) 


we can make the following explicit choice of the polarization basis, 




V2 


(0, cos 9 cos <p + i sin <p, cos 9 sin <p ± i cos <p, — sin 9 ) . 


(19) 


Therefore we choose the polarization vectors to be helicity eigenvectors. For the special choice q M = (po, — p) 
they satisfy (fTHl) and (fT71) . 

We denote by J*a 3 ,a 4 ;Ai,a 2 the helicity amplitudes for the scattering process 1 + 2 -+ 3 + 4 of spin 2 
particles with helicities A $. The helicity amplitudes satisfy a number of relations as a consequence of parity 
conservation and time reversal for a reaction a + b — > a + b, Bose symmetry for a reaction a + a —> b + c, 
and invariance under particle-antiparticle conjugation for a reaction a + a — > b + b, namely 

7*A 3 ,A 4 ;Ai,A 2 ( t) A 3 , —A 4 ; —Ai, —A 2 5 

-^A 3 ,A 4 ;Ai,A 2 ( f) -^Ai,A 2 ;A 3 ,A 4 ? 

-^A 3 ,A 4 ;Ai,A 2 ( f) ^A 3 ,A 4 ;A 2 ,Ai 5 

-Fa 3 ,A 4 ;Ai,A 2 = (-1) A M 7*A 3 ,A 4 ;A 2 ,Ai , (20) 


where A = Ai — A 2 and p = A 3 — A 4 . 

The most significant advantage of working with such amplitudes is that they are defined as Lorentz 
invariant quantities for massless particles, which implies we are automatically computing gauge invariant 
quantities, only dependent on the Mandelstam invariants s = — (pi + P2) 2 , t = — (pi — P3) 2 
and u = — (pi — P4) 2 . 

Using relationships G 0 , out of 16 amplitudes in four dimensions, we can single out four independent 
helicity amplitudes [5| 


A ( + + , ++) = i*2,2;2,2, A (H , -I ) = F 2 - 2;2,-2, 

A (++, H ) = 7*2,2;2,-2, A (++,-) = 7*2,2;—2,-2- 
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2.2. The two-point function 

Let us start by considering the propagator of the Einstein-Hilbert theory assuming the expansion of the 
metric around a flat background 


9fiv — ‘hn 1 ' 4 " d P 


( 22 ) 


where h pv is the fluctuation field. We will denote by the n-th order term of the expansion of operator O 
in the fluctuation field. Furthermore we will distinguish between on-shell and off-shell fields by underlining 
the latter. 

The equation of motion is given by 


Seh = — ^ k d d x y/—gR 


^EH 

S 9, 




Sg P , = -2 k~ d 2 fd D x V^g {^f VR ~ Sg„ 


(23) 


We can always assume ((SSeh/^AVu/)^ = 0 because R= R P J = R^ pa = 0 as a consequence of having 
chosen a flat background metric. The two point function is obtained by considering the expansion to the 
second order 


<A 2 ) — _Pk- 
j eh — 


-2 


' d D x(V=g R ) {2) = -2k 


-2 


d°i 


(v^) (2) RW + (v^) (1} R (1) + {V=g) m R (2) ■ (24) 


,(!) 


( 0 ) 


For an on-shell graviton 


R w = d a d b h ab -Dh a a = 0 , 

R$ = l (-d P d„h a a + d p d a h av + d„d a h ap - □ V) = 0, 


R$eo- = \ {-d^dpKo- + dpdaKp + d v d p h pa - d v d a hp. p ) ± 0 , 
which in particular imply the minimum condition for the quadratized Einstein-Hilbert action 
<5 2 Seh 


\ (o) 


0. 


(25) 


(26) 


,Sg a p(x)Sgj S (y) J 

For the higher derivative theory (1131) we need to add the second order expansion of the quadratic terms 
Squad 


Squad = [d D x (^) (0) ( 7o i? (1)2 + 72^ )2 + 74^) • 


(27) 


To be able to define the inverse of the kinetic operator we add to the Lagrangian the standard gauge-fixing 
term 


S gf = -2 ^ Jd D x {d a h av - /3cU“) 2 , 


where we choose a = 1 and (3 = 1/2, corresponding to the usual harmonic gauge condition 

d a h av = d v h% . (28) 

This is consistent with the physical conditions (fTTl) chosen for on-shell gravitons. We can rewrite the kinetic 
operator in the momentum space as follows, 

\^- k) (s£&y Sy l0 ^ (t) 55 ■ < 2! » 

A standard procedure can be employed to compute the propagator iO~g ^g{k) as reviewed in | Appendix A 










2.3. The three and four-graviton vertices 

Three graviton vertex — The three graviton vertex is obtained by collecting all terms in the action that 
are cubic in the graviton field h pl ,. Since we are interested in the on-shell four-point function, in particular 
in the scattering amplitudes of four massless gravitons, we do not need to go through a lengthy computation 
keeping track of all terms eventually vanishing on-shell. Therefore, we can make full use of the simplification 
brought about by the linearized vacuum equation of motion for the physical field h in the harmonic gauge 
(1281) . i.e. Oh^ = 0. Actually, our choice of polarizations is such that we can assume the conditions m 
all along our computation, which greatly simplifies the algebra. In fact, these conditions imply that all the 
scalar operators are vanishing on-shell at linear order in h^, including the scalar curvature i?) 1 ) and the 
root of metric determinant y/—g 1 ' 1 ' 1 ■ One can further show that R$ = 0 due to the linearized EOM. The 
diagrams we need to compute are such that for each three graviton vertex two gravitons out of three are 
on-shell and one is off-shell because it must be contracted with the propagator of the internal line of the 
diagram. This means we can consider the physical and off-shell gravitons as different fields, even from the 
combinatorial point of view. Hence the procedure is first to expand the action to the first order in the 
off-shell field h jJV and then to the second order in the physical field h pl/ . 

It is convenient to rewrite the action CE3D in terms of the only two combinations of the 7 £ parameters 
that appear in the propagator, namely 

S g = —2 k d 2 Jd D x y/^g^R + (70 — 74 ) R~ + (72 + 474 ) R 2 V + 74 (-ffjLpo- — 4 R^ v + R ~) ^ . (30) 


We note the last term is the famous Gauss-Bonnet density, which is topological in four dimensions whereas 
for generic higher dimensions it gives rise to vertices only. 

For the Einstein-Hilbert action Seh = —2k)) 2 IdP Xy/ —gR we simply have for the three graviton vertex: 


./^ eh \ {2) . 


— * ( 2k d 2 ) 




V %«/ J 

Similar expressions arise from the scalar curvature square action 
S’o = - 2 k )) 2 (70 - 74 ) Jd° x yfPjR 2 , 
the Ricci square action 

S' 2 = -2k)) 2 (72 + 474) J d D x y/^R] 


1 ( 2 ) 

(^g^R-R^A = i (-2k), 2 ) (lgf v R { ?) - R^ 2 ^ (31) 




and the Gauss-Bonnet action 


S' A = -2k~ 2 74 J d D x yf^g (R 2 wpa - AR% + R 2 ) . 


The outcome is: 

\( 2 ) 

y/w 


"I (2) 


2v /3 5 ( /"V 2 - + - g^R - R! u/ ) R 


( 2 ) 


(32) 


(33) 


(34) 


‘ (£)'=* ( - 2 ^ 2 ) ( 7 .-„) 

= 2 i (-2 k~ d 2 ) (70 - 74) (»r □ - ^d v ) RW , 

i = * (-2 k- d 2 ) (72 + 474 ) v=g (\g^R KX R K x - 2 R^R\^ 

+* (—2k)) 2 ) (72 + 4 74 ) [yf^ (-g^g^V 2 - ^"V K V A + g^W x V" + g vX V K W) R kX 
= i (-2k)) 2 ) ( 72 + 474) (-7 fV A D - ri^d K d x + ■ n ^d x d v + 7 f x d K d») R (2 l 


(35) 


( 2 ) 


(36) 
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where □ is the flat d’Alambertian operator. In D = 4, as far as the computation of vertices is concerned, 
we can ignore the presence of £>4 because of the Gauss-Bonnet theorem. The terms in (I3T1) . (1351) and (1331) 
contributing to the three graviton vertex with two legs on-shell are expressed in terms of the following second 
order expansions 

R {2) = -d b h ac d c h ab + 1d c h ab d c h ab , 

R$ = \d^,h ab d v h ab + h ab (d b d a h^ + d^d v h ab - d b d^h ua - d b d„h^ a ) + d b h^ a ( d b h va - d a h ub ) . 

Four-graviton vertex — For the on-shell four-graviton amplitudes we only need the four-graviton vertex 
with all the gravitons on-shcll. Therefore, by the same arguments used for the three vertex, we can argue 
that 

* (W 4) = -2 Jd D x {V=gR) {4) 

= -2 Jd D x (i? (4) + {V=d) {2) R {2) ) , 

i ( S'o) {4:) = -2 Kd 2 (70 - 74) Jd D x , 

i (S') (4) = -2 k~ 2 (72 + 474 ) J d D x (r$) 2 . (37) 

Apart from (y/—g)^ = the only new quantity appearing in the four vertex is R^\ whose 

expression is given in | Appendix B| We remind the derivation of all vertices from S'eh, S' 0 and S' 2 is valid in 
any dimension D. 

3. The amplitudes in Stelle gravity 

In this section we explicitly evaluate the four-graviton scattering amplitudes for quadratic Stelle gravity 

using the method of Feynman diagrams. We distinguish the two cases of four and higher dimensions. In 

the first case we have already introduced all the ingredients which will be used in the computation. For the 
second one we will add new vertices and extend the form of polarization tensors. This second case will be 
the first extension of our results first obtained in quadratic four-dimensional Stelle theory. 

3.1. J^D Stelle gravity 

In Born approximation there are four diagrams to be considered, which include the contact one and the 
ones with virtual propagation in the s, t and u channels. They are depicted in Fig[lJ while the kinematics 
of the process is showed in Figj2] 

The exchange diagrams can be computed by contracting two three graviton vertices with the propagator. 
Out of the five possible tensor structures in the propagator (see | Appendix A| ) only three (Xi, X 2 and X 3 ) 
are seen to have a non-vanishing contribution as a consequence of the fact that we are considering on-shell 
gauge invariant amplitudes. In order to be able to carry out the algebra effectively, it’s very convenient to 
express the Mandelstam variables in terms of the energy E and the scattering angle 9 in the center-of-mass 
reference frame, 

s = AE 2 , t = -2E 2 (1 - cos 9) and u = -2E 2 (1 + cos 9) . (38) 

With this choice of variables and assuming all the momenta lying in the plane defined by <f> = 0, the incoming 
(pi and p 2 ) and outgoing ( p 3 and p±) gravitons are identified by the momenta (1181) and polarizations (1151) 
obtained from (ED, where 9 is chosen to be 0, 7 r, 9 and 7 t + 9 respectively. For outgoing gravitons the 
complex conjugates of polarization tensors should be used. 


10 






s channel 


t channel 


u channel 



contact 


Figure 1: Tree-level Feynman diagrams giving rise to the process of four-graviton scattering. 


The contributions from different diagrams to the amplitude for gravitons with all helicities +2 are 

A s (++, ++) = -i (- 2^4 2 ) ^E 2 sin 2 0 - i (- 2 k 4 2 ) ^-E a 8(70 - 74 ) + (72 + 474 ) (cos 20 + 3) 

o 4 L 


At (++, ++) — — 


. (— 2k 4 2 ) cos' 


,4 6 
2 


E 2 (392 cos 0 + 80 cos 20-8 cos 30 - 720) 


A u (++> ++) — — 


256(cos0 — 1) 

+(70 — 74 )E 4 (4cos0 + 16 cos 20 — 6 cos 30 — 4 cos 40 + 2 cos 50 — 12) 

+(72 + 474 ) 1? 4 (2146 cos 0 — 216 cos 20 — 99 cos 30 + 6 cos 40 + cos 50 — 1838) 
. (— 2k 4 2 ) sin 4 | 


E 2 (392 cos 0-80 cos 20-8 cos 30 + 720) 


Ac 


(++, ++) — i 


256(cos0 + 1) 

+( 7 o — 74 )-E 4 (4 cos 0 — 16 cos 20 — 6 cos 30 + 4 cos 40 + 2 cos 50 + 12) 

+(72 + 4 7 4)£" 4 (2146 cos 0 + 216 cos 20 - 99 cos 30-6 cos 40 + cos 50 + 1838) 

(-2. 4 - 2 ) r 


E 2 (160 cos 20-8 cos 40 + 872) 


(39) 

(40) 


(41) 


(42) 


1024 

+(7o — 74 )E i (—34 cos 20 + 4 cos 40 + 2 cos 60 + 18460) 

+(72 + 474)£ 4 (1711 cos 20 - 46 cos 40 + cos 60 + 9598) 

These results themselves contain a nontrivial piece of information. In fact, in such a covariant gauge as the 
harmonic one we have chosen, at any energy scale the propagator behaves like 


1 


1 


1 


k 2 (1 — 7 'k 2 ) k 2 k 2 — 7 ,_1 


(43) 


and because of the small k expansion [ k 2 (l — 7 'fc 2 )] 1 = k 2 + 7 ' — 7 12 k 2 + ... we would expect quantities 
dependent on arbitrary powers of E 2 , whereas we find that only energies up to the fourth power show up 
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Figure 2: Kinematics of the process hh —> hh. 


and in particular no energy powers appear in the denominators. This is the sign of nontrivial cancellations 
happening between the three graviton vertices and propagators for each diagram. The ghost poles in (1431) 
disappear at any energy scale and moreover we are safe from IR divergences. However the contribution from 
each diagram is not separately gauge invariant and only the sum giving the full amplitude matters. 

Actually, it’s easy to check that the full scattering amplitude reads 


A ( + + , ++) — As (H—b, ++) + At (+ + , ++) + Au (++) ++) + -^contact (++> H—I") 



which is the same result as the one for the tree-level amplitude of Einstein theory. We can also find a similar 
result for the other three independent helicity amplitudes 


A (-1—, 4—) 
A (++, H—) 




(1 + COS 0) 4 

sin 2 9 


A (4—(-, —) — 0. 


(44) 

(45) 


As a consequence of another nontrivial cancellations these results coincide with the very well known one Q 
for tree-level graviton scattering in Einstein theory in spite of the fact that the dimensional arguments that 
completely constrain the form of helicity amplitudes in the latter case [fj can no longer be naively applied 
to quadratic gravity in D = 4, because of the presence of new dimensionless couplings in the theory. We 
will find a natural explanation for these surprising results in section [5j 


3.2. D > 4 Stelle gravity 

For D > 4 the situation becomes more complicated. First of all, the Gauss-Bonnet term is no longer 
topological and so it gives a non-vanishing contribution to the three and four-graviton vertices. Apart from 
the terms we have already computed, for S 4 = (—2k2) 2 ) 74 Jd D x^J—gR? pvpG we find 


(|p) ( = (—2 k ~ d 2 ) 74 2^ (jg^R KXp,7 R K x P * - g UT R pXprJ RrXpa + VaV^ 1 * 1 ^) 


( 2 ) 


-2^ 2 ) 74 [fT (Xi.) - 2 ri VT ri^R% a R { X, f^R™ + UR^\ 


( 2 ) 

\vpG 


i (S 4 ) (4) = (- 2 ^) 74 J d D x[ (R ^) 2 + 2 R$ pa R$ pa +8 g^R^R 

+ i2^(i) 5 -(i) i? (i) p<Ti? w a + 4 ^ 2 ) 41 ^^ + v^ (2) 


2 1 


(46) 
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The three and four-graviton vertices are therefore determined in terms of the on-shell non-vanishing quanti¬ 
ties, whose explicit expressions are given in |Appendix B| As we will see from the results below, it is crucial to 
include the term — gR 2 ^ pa in the action in higher dimensions. This term is naturally motivated by effective 
field theory considerations and addition of it is in another direction, how we can extend four-dimensional 
theory. 

An additional complication is given by the fact that in D > 4 the graviton has D( ' L> 2 ^ polarizations. 
This is a trivial consequence of the fact that the little group for a massless particle in D dimensions is 
SO(D — 2) and that the graviton is by definition identified with the traceless symmetric rank 2 tensor 
representation. Because a massless vector particle with momentum ( p 2 = 0) has (D — 2) independent 
polarizations ej/(p, A) we can choose the Lorentz gauge condition p^e^ = 0 and identify up to a gauge 
transformation + 7 (for any 7 ) to single out an irreducible representation of SO(D — 2) with A 

identifying the elements of the vector basis. This basis can be conveniently chosen as the one such that 

JPe>,A) = A*e>,A), 

with H l the elements of the Cartan subalgebra , i = 1 ,...,r. For SO(D — 2) we can break up the 
(D — 2)-dimensional space into [D/2] — 1 different two-dimensional subspaces. The rotation operator W k 
(k = 1,..., [D/2] — 1) acts on the subspace (2k, 2k — 1). Therefore polarization vectors are identified by the 
D — 2 weights of the fundamental representation of SO(D — 2). If we assume the momentum p M to be in 
the (D — l)-th spatial direction and W^ D D\- X to be the generator of rotations in the (D — 3, D — 2) plane, 
then we should consider in particular the polarizations corresponding to the [D/2] — 1 dimensional weights 
(0,..., 0, ±1), i.e. the lA-vector (0,..., 0,1, ±i, 0). In the following computation we will concentrate on this 
structure of D-dimensional polarization vectors. Now we want to repeat the discussion in the sectioned] for 
a general dimension D. 

For a generic momentum whose spatial components are only in ( D — 3, D — 2, D — l)-subspace 

= (.Po, 0, ..., 0, po sin 9 cos (/>, po sin 9 sin </>, po cos 9 ) , (47) 

we have 


(p, ±) = (0,..., 0, cos 9 cos (j) =F * sin </>, cos 9 sin (/ ± i cos </>, — sin 0) . 

These polarizations can be used to construct two traceless symmetric tensors satisfying the gauge condition 
P^e^v = 0 and equivalent up to a gauge transformation e^ —> e^ + a^p u + a^p^, where a ■ p = 0, 

±) = (P> ±) e ?0. ±) ■ 


Gravitons in D dimensions, which we consider here, are with two possible polarizations + or —, similarly to 
the case in four dimensions. 

In this framework we have carried out the computation for the amplitude A D (++, ++) in D = 5 and 6 
dimensions. These are the results 


A D ~ 5 (++,++) = 

A D = 6 (++,++) = -i\ 


16 F; 6 7 | [1 + 8 E 2 ( 3(70 - 74) + (72 + 4 y 4 ))] _ 2 1 \ 

(1 - 4£ 2 (72 + 474)) [3 + 4 E 2 ( 16(70 - 74) + 5(72 + 474))] sin 2 9 j ' 

8F; 6 7| [1 + 8 E 2 (3(70 - 74) + (72 + 474))] _ 2 1 \ 

(1 - 4£’ 2 (72 + 474)) [1 + 2 E 2 (10(7o - 74) + 3(72 +474))] sin 2 9 j 


Quite beautifully, each amplitude is the sum of two terms, one comes from the usual Einstein gravity, 
whereas the other one is overall proportional to 74 . This implies that in the absence of the Ricmann square 
action S 4 the (quite boring) result found in four dimensions holds true in higher dimensions too, making 
the scalar curvature square term So and the Ricci square term S 2 undetectable in the tree-level graviton 
scattering amplitudes. This is why the presence of the S 4 term in the action is crucial in higher dimensions. 
The mechanism of cancellation in the former case is the same as in four dimensions. 


13 









Another interesting feature is that, apart from the standard Einstein term, the remaining dependence 
of the amplitudes on the parameters of the quadratic terms is only through the combinations 70 — 74 and 
72 + 474. To better understand we first define a new quadratic action for gravity with parameters 7 q, 72 
and 74 according to 


S g — 2 k d 



+ 7o-^“ + 72-ftjL + 74GB^ , 


( 50 ) 


which is completely equivalent to (USD, and where GB denotes the Gauss-Bonnet Lagrangian. We find 
precisely that 7 ' 0 and 72 are equal to combinations appearing in (l 48 l) and (l 49 l) . Therefore those results can 
be rewritten in a more compact form 


a d -5 ( ++ ++) = -i2- 

"5 


16 E 6 7 | [l + 8 E i2 ( 3 7 q + 72)] _ 2 1 \ 

(1_4£ ,2 7 ')[3 + 4^; 2 (167'+57')] " sin 2 0j’ 


A D ~ 6 (++,++) = -i\ 
«6 


8 -E ’ 6 7 4 [l + 8E 2 (3^q + 72)] 2 1 1 

(1 - 4 £2 7 ') [1 + 2^(107' + 37')] sin 2 9 j ’ 


( 51 ) 

( 52 ) 


The dependence on the parameters 7$ in the two above formulas can be easily explained diagrammatically. 
We know that the propagator derived from ( 15 ( 71 ) doesn’t depend on 74 coefficient in any D, only cubic and 
quartic vertices derived from the Gauss-Bonnet term possess such dependence (precisely expressions for them 
are linear in the 74 parameter). The propagator depends only on k^ 2 , Y 0 and 72. Hence our conclusion is 
that besides contribution to amplitudes m and ( 1521 ) from Einstein gravity, we have additional contributions 
from exchange diagrams, where both three graviton vertices are derived from the Gauss-Bonnet term and 
on the internal line we have the full propagator of the theory. The dependence on 7^ and 72 is only through 
the propagator, but not through vertices. The dependence on k^ 2 is determined by dimensional reasons. 
The fact, that there is no any linear dependence on 74 in the final results forces us to believe that here 
we are witnessing another cancellation between contact diagram with vertex from Gauss-Bonnet term and 
exchange diagrams with two different vertices (one from GB, the second one from standard terms R , R 2 
or i? 2 „). This interpretation of the results is quite natural, because we have modified the theory (with 
already many cancellations) only by addition of new vertices coming from the Gauss-Bonnet term, but the 
propagator has remained the same. We have also computed the amplitudes for other choices of polarizations 
and checked that their general properties are similar. 

Moreover, as it is obvious from ED and ( 1521 ) the new terms in the amplitudes are associated with the 
appearance of arbitrary powers of E 2 in the expansion of the denominators in the infrared regime, while 
in the ultraviolet regime the highest power is of course E 4 , because this is a quadratic gravity. This result 
already extends the previous findings [ 3711571 ) , which were restricted only to the order E 4 in energy expansion 
around E = 0 . When we treat the theory m as fundamental, then we do not need to focus on the low 
energy limit of the amplitudes and we have the exact energy dependence in the tree-level amplitudes. By 
comparing the cross sections computed in [ 57 ], we could in principle read out the values of the parameters 
7p, 72, and 74. Then it would be natural to associate the parameter 74 with the strength of Gauss-Bonnet 
interactions, while two other "f'o and 72 would be related to the masses of the ghost and the curvaton, which 
appear in the spectrum of quadratic gravity. However, this theory is not unitary (due to the presence of the 
ghost) and we cannot conclude about new physically meaningful contributions to the graviton scattering 
compared to the amplitudes computed in Einstein gravity. In the next section we would like to address the 
same issue in a more realistic theory. 


4. Four-graviton scattering amplitudes in nonlocal gravity 

In this section we explicitly calculate the four-graviton scattering amplitudes for the weakly nonlocal 
theory © with zero potential (V = 0 ), namely 

£ g = - 2 k~ d 2 ^[R + R 10 (U)R + R^ 2 (D)R^ + R^ pa . ( 53 ) 
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This is another direction of extension of our original four-dimensional result in Stclle gravity from 13.ll Here 
we assume the theory to be valid in any dimension D, but the explicit formulas will be given for the case 
D = 4 (in higher dimensions only the numerical coefficients change, but the final results are the same). 
Additionally results from this section can be used in particular for any local higher derivative theory, whose 
action is quadratic in gravitational curvature and contains a finite number of derivatives. We remark that 
for below results to hold true the assumptions about super-renormalizability (or only renormalizability), 
unitarity and nonlocality of the theory are not by any means necessary. Since we compute the tree-level 
scattering these issues are irrelevant. We decide to speak about the theory ©, because this is a candidate 
theory to have a good behavior at quantum level too and could be viewed as a fully consistent realization 
of quantum gravity in quantum field theory framework. 

The stunning result is once again the same of Einstein gravity for the case 74 (D) = 0 in any dimension. 
Since on-shell the expansion of two the simplest curvature invariants is R ~ 0(h 2 ) and Ric ~ 0(h 2 ), the 
form factors are spectators in the expansion in number of gravitons and many of the results of the previous 
section still apply to the general nonlocal theory (1531) . In more technical terms, when we compute the 
variations to find propagator and vertices, we don’t need to vary the form factors and original covariant 
boxes □ being the argument thereof. Putting all the n-point functions on flat spacetime we can substitute 
these arguments with flat d’Alembertian operators. Finally after going to momentum space we can easily 
replace them with negative invariant squares of the momenta - Mandelstam variables — s, t and u for each 
channel respectively [18|]. Only for the contact diagrams we have to be careful and we closely investigate on 
which graviton legs the operator inside the form factor acts to properly associate dependence on s, t and u 
in form factors in momentum space. 

We here report the result in D = 4 and with form factors appearing only inside curvature scalar and Ricci 
tensor squared terms (no terms with Riemann tensors sandwiching the form factor 74(D)) for the amplitude 
M(+-b, ++) and we explicitly show the cancellations leading to the important result announced above. The 
amplitude A(+- b, ++) gets contributions from the contact diagram and ones with graviton exchanges in the 
s, t and u channels (compare Fig[T]), namely 


A s (-|—b, H—b) — —2ft 
At (-1—b, H—b) = —2ft 4 2 I — 


_2 / 9 t(s + t) 9 


+ t^7 2 0) (s 2 + (s + 2 1 ) 2 ) + ^s 2 7 o(s)^ 


4 ' 8 s 

1 (s 3 — 5 s 2 t — st 2 + t 3 ) (s + t) 
8 sH 


1 (2s 4 — 10s 3 f — s 2 t 2 + 4sf 3 + f 4 ) (s + t) 2 1 t 2 (s + t) 


+ 7772W 

16 5 - 

o / 1 (s 3 — 5 s 2 u — su 2 + u 3 ) (s + u) 2 

Au(++,++) = - 2ftJ 2 ~-- 


o7o(i)- 


c 4 > 


(54) 


(55) 


1 (2s 4 — 10s 3 u — s 2 u 2 + Asu 3 + u 4 ) (s -b u) 2 1 N u 2 (s-bu) 4 

+7772(u)^--4- L -- + ^7o(») o4 


16 


_ 2 1 1 s 4 + s 3 t — 2 st 3 — f 4 9 , N / o . , 9 


- 4 C ontact(+-b, + + ) — — 2ft 4 ( — - ^ 


-7^72(s) (s 2 + (s + 2f) 2 ) —-s 2 7o(s) 


1 (2s 4 — 10s 3 t — s 2 t 2 + 4sf 3 + t 4 ) (s + t) 2 1 t 2 (s + t) 4 

1 . (2s 4 — 10s 3 u — s 2 u 2 + 4su 3 + u 4 ) (s + u) 2 1 ,« 2 (s + «) 4 ' 

~16 72(U) -?-g 7o (u) ^4 


(56) 


(57) 


The full amplitude is given by the sum of above contributions (1541) -b (1551) + (1551) + (1571) . and the result is: 


-A( + + , ++) — M s (+ + , ++) + At{- b+) + + ) + Au{- 1 —b, ++) + 7 lcontact(+ + ) ++) — - 4 (-|—b, ++)eH 


, ( 58 ) 
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where „4(++, ++)eh is the amplitude for the Einstein-Hilbert theory. Notice that all terms but the first of 
equation (1571) cancel with the last two terms of the first line of (1571) . Analogously terms from (1551) cancel 
with the second line of ED and those from m cancel with the last line of m- For the other helicity 
amplitudes (ED we get analogous simplifications. 

Once the operator quadratic in Riemann tensor (or Weyl tensor, or the generalized Gauss-Bonnet term) 
is turned on the scattering amplitudes change radically. Not only amplitudes will explicitly depend on the 
form factor 74(D), but also on the other form factors in the theory, namely 70(D) and 72(D), as it is evident 
from the five- and six-dimensional results reported in ( 1511 ) . ( 1521 ) for the case of constant form factors. There 
74 is an overall factor, and if it does not vanish, the results will depend also on 70 and 72. However the 
mechanism for acquiring such dependences is exactly the same as described in section 13.21 The case of 
computation with non-vanishing form factor 74(D) is qualitatively different from presented here. In this 
situation to obtain vertices we need to vary the form factor and covariant boxes in it. The reason for this is 
simple to explain the expansion of Riemann tensor on-shcll starts at the first order, Riem ~ O(h), which 
is in opposition to the case for Ricci tensor and curvature scalar. 

When the nonlocal generalized Gauss-Bonnet operator 

GB 74 ( D ) = Riem74(D)Riem — 4Ric74(D)Ric + Ry4(D)R ( 59 ) 

is switched on the amplitudes will depend explicitly on 74(D) (through vertices) and on 7q(D) = 70(D) — 
74(D) and 72(D) = 72(D) + 474(D) (through propagators). This operator is non-trivial in any dimension 
D > 4 , and gives rise only to new vertices, but not to the full propagator of the theory. As it is now 
obvious from the discussion in section 1X71 the vertices derived from this term contain derivatives of the form 
factor with respect to its argument (up to the second order). Moreover the argument of the form factor, 
which is the covariant box operator (acting on a tensor field up to rank four) must be varied too (up to the 
second order), which complicates the situation quite a lot. Therefore the computation of such vertices is 
quite involved and we don’t attempt to present the exact results here. We only remark that the results for 
constant form factors in section I5~21 for D > 4 are consistent with these qualitatively described here. The 
presence of a non-constant form factor 74(D) in a nonlocal theory ( 1531 ) is crucial even in D = 4 and there in 
an ideal graviton scattering experiment we could fit the form factors 70(D), 72(D) and 74(D) by measuring 
the cross section. 


5. General n —graviton scattering amplitudes in local and nonlocal theories 


Having discussed the scattering amplitudes in four-dimensional Stelle theory and in simple extensions 
of it in the previous sections (higher dimensional setup, inclusion of terms quadratic in Riemann tensor, 
and nonlocal form factors) now we wish to explore other possible directions about how we could extend our 
results. The main motivation is to seek for a well-defined setup in which amplitudes will differ significantly 
from those obtained in Einstein theory and therefore would permit for unambiguous verification of the theory 
by comparing predictions for cross sections with hypothetical experiments on gravitational wave scattering. 
As we will see in this section it is still difficult to depart from very ubiquitous results of Einstein gravity. As 
a first step we will try to describe the situation when n-point correlation functions are considered. In this 


area the theorem (first proven by Anselmi in [32l . 1331 ]) will reveal to be very enlightening for the tree-level 


situation. We will again find typically only standard results from Einstein theory and a justification for this 
will be given. A special role will be assigned to the Riemann tensor and to the form factor 74 (D). At the 
end we will comment on other possible extensions by the inclusion of other operators: local (and higher in 
curvatures) and in a sense more nonlocal, and by going beyond tree-level. In this way we will touch on any 
reasonable extension of the original theory m- Exploration of all these directions in an exhaustive way is 
beyond the scope of this paper and we will leave it for future publications. 

Let us consider the case of higher than four-graviton scattering amplitudes at tree level. For this goal 
we present here the following theorem. 

Theorem. All the n-point functions in any gravitational theory (in particular super-renormalizable or finite) 
with an action 


-2 


jCgr — 2k j , 


\f--g [R + Ryo(D)R + Ric 7 2 (D)Ric + V(R, Ric, Riem, V) ], 


(60) 
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give the same on-shell tree-level amplitudes as the Einstein-Hilbert theory, £eh = —2 nf/ provided 

that the potential V is at least quadratic in Ric and/or R. In particular for any theory in which we can 
recast the potential in the following form 

V = Ric V Ric = R p /V{ R, Ric, Riem, V)] pvpa R pa , (61) 

the theorem is still valid (V is in full generality a differential operator with contravariant indices /i, v, p, a 
acting on the Ricci tensor to the right, containing at least one power of gravitational curvature.) 

Proof. The proof is based on the field redefinition theorem proved by Anselmi [ 32 } at perturbative level 
and to all orders in the Taylor expansion of the redefinition of the metric held. 

First we assume that we have given two general weakly nonlocal action functionals S'(g) and S(g'), 
respectively defined in terms of the metric fields g and g', such that 

S’( 9 ) = S(g) + E i (g)F ij (g)E j (g) , (62) 

where F can contain derivative operators and Ei = SS/Sgi is the EOM of the theory with action The 

statement of the theorem is that there exists a held redefinition 

g'i = gi + AijEj A ij = A ji, (63) 

such that, perturbatively in F , but to all orders in powers of F, we have the equivalence 

S/g) = S(g'). (64) 

Above A ^ is a possibly nonlocal operator acting linearly on the EOM Ej, with indices i and j in the held 

space, and it is dehned perturbatively in powers of the operator Fi/g), namely A l:j = Fij(g) + .. . Let us 

consider the hrst order in the Taylor expansion for the functional S(g'), which reads 

S(g') = S(g + A g) k, S(g)+ j^Sgi = S(g) + E l 6g i . (65) 

If we can hnd a weakly nonlocal expression for Sgi such that 

S’(g) = S(g) + Ei Sgi ( 66 ) 

(note that the argument of the functionals S' and S is now the same), then there exists a held redehnition 
g —> g' satisfying HMD- Hence the two actions S'{g) and S(g') are tree-level equivalent. □ 

As it is obvious from above, in the proof of our theorem it was crucial to use classical EOM Ei. In 
the theory (1601) this implies Ric = 0 in vacuum regions without the presence of any matter source. These 
are the conditions we have imposed on the linearized level when we dehned asymptotic states of on-shell 
gravitons as perturbative states of our theory. 

Now we can apply the above held redehnition theorem to our class of theories ©> where we do not 
include terms with Riemann tensor Riem. Since we are interested in S(g') = SehG/) and S'(g) = S gI (g), 
the relation (1521) reads 

S(g') = Seh(3) + RMF^ p %g)R pa {g) = S\g). (67) 

Here we also used that in the spectrum of Einstein-Hilbert theory we only have the massless spin 2 graviton 
(contrary to the case of the theory described by the action (1501) with polynomial form factors) and the 
on-shell scattering of such particles we relate in the two theories. It would be clearly nonsensical to apply 
the theorem for scattering of other particles (appearing e.g. in theory (1551) 1 and attempt to relate it to the 
scattering amplitudes in the Einstein two-derivative theory. 


2 Here we use a compact deWitt notation and with the indices i, j on fields we encode all Lorentz, group indices, and the 
spacetime dependence of the fields. Additionally, we assume that the field space is flat and we do not need to raise indices in 
sums there. 
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If the potential V is at least quadratic in R and/or Ric, namely takes the form (1611) . where V is a rank 
four tensor made of any tensor including Riemann tensor, its contractions and derivative operators, then 
the tensor Fij(g) used in the field redefinition (1621) exists, it is weakly nonlocal, and equals to 

F pv ' pcT = g^g p,T 70 (D) + gpp g va 72 (D) + V(R, Ric, Riem, Vf vpa . ( 68 ) 


In the additional part of our theorem we will discuss about the possibility of recasting the potential V in 
the form m- Of course due to the known ambiguity related to an order of writing covariant derivatives 
(non-commuting on a general manifold) in tensorial expressions, this last remark is not very precise. More 
precisely we require that the contribution to the EOM from the potential V should vanish, when the ansatz 
Ric = 0 is used. (This is the conclusion about vacuum spacetime in Einstein gravity, which on the linearized 
level coincides with the on-shell conditions for the massless gravitons.) We do not say that in order to use 
the theorem the initial potential V must be in the form (1611) . it only must be possible to cast it in such 
form. 

Our results about Stelle theory in four dimensions reported in section 13.11 can be understood in the 
following way. First, we may employ Gauss-Bonnet theorem to reduce the action to the form with only 
the Ricci scalar and the Ricci tensor square terms. Secondly, we can use the field redefinition theorem to 
prove that this theory is equivalent to the Einstein-Hilbert theory regarding the on-shell tree-level graviton 
amplitudes. Another convenient choice is to start with a theory written in a Weyl basis with the quadratic 
part consisting of the Ricci scalar square R 2 and the Weyl tensor square C 2 . We can take the limit in 
which only the coefficient in front of C 2 survives so that the theory is now conformally invariant at classical 
level for D = 4. The operator C 2 should contribute with the fourth power in the energy to the four- 
graviton scattering amplitudes, but conformal invariance requires the scattering amplitudes to be numbers 
independent on the scale. Therefore, the amplitudes must be zero because the graviton field is dimensionless 
and there is no other scale in Weyl gravity. On the other hand, in J\f = 4 super-Yang-Mills theory we can 
have non-zero amplitudes because the gauge bosons have energy dimension one. However, we explicitly 
proved that also the operator R 2 does not give any contribution to the amplitudes, so that we can conclude 
the scattering amplitudes for any purely quadratic gravity in D = 4 are vanishing. 

As a special case, our explicit computation confirms that the four-graviton scattering amplitudes in four 
dimensional Weyl conformal gravity [52| is identically zero, namely 

C = —a g y/^g C^ ptT C p,l ' pa => M(4-graviton) = 0. (69) 


This derives from the fact that in any CFT S matrix must be trivial. 

In D > 4, operators quadratic in the Riemann tensor (with or without form factors) cannot be recast in 
the form of operators quadratic in the Ricci and scalar curvatures without introducing extra vertices because 
the generalized Gauss-Bonnet operator (1591) is no more topological. Therefore, they could contribute to the 
scattering amplitudes. This completely explains the results found in section f3.2l and also explain why in the 
case of form factors sandwiched among Riemann tensors the results will deviate from the ones computed in 
Einstein theory. In the four dimensional case the generalized Gauss-Bonnet term (1591) does have an impact 
on the EOM and the Ricci-flatness ansatz is not valid in the vacuum of the theory. Therefore, also in this case 
we cannot apply the theorem and the amplitudes differ from those of the Einstein theory in agreement with 
the results reported in section 01 In particular the amplitudes in the Einstein-Hilbert action supplemented 
by the two-loop correction term (Goroff-Sagnotti term [36j), understood as effective field theory, will not 
coincide with those from pure Einstein theory. 

We comment also on applications of the above theorem to the case of dimension D. For the case of a finite 
theory the result can be extended to any order in the loop expansion if we neglect the finite contributions to 
the quantum action. In even dimension we can easily achieve finiteness by using of killers constructed at least 
out of two Ricci tensors (compare formula © for D = 4). Note that in odd dimension, for 7 > (D — l)/2 
in the description after formula the theory (15H1) is finite in DIMREG without the need to add any killer 
operator. We also note that in D = 3 the Riemann tensor is not independent from the Ricci tensor and the 
scalar, therefore three-dimensional Einstein-Hilbert gravity without matter is finite at quantum level (but 
without perturbative degrees of freedom). 
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For the case of a 1-loop super-renormalizable theory in D = 4 the theorem can be applied at any 
order in the loop expansion including quantum loop divergences of amplitudes and the quantum logarithmic 
corrections coming together with the one-loop running of gravitational couplings. However, we again expect 
deviations from Einstein Hilbert amplitudes due to other quantum finite contributions. It is crucial here that 
the theory is one-loop super-renormalizable, because we only have divergences at the controllable one-loop 
level. In a general renormalizable theory we would have divergences of the type R 2 and R 2 pv at any loop 
order in D = 4, and the structure of the RG equations for the running coupling constants or the finite 
terms in the effective action will be much more complicated. In the case of theories renormalizable and 
(super-renormalizable) in higher dimension D > 6 we cannot apply our theorem any more, because then 
operators of the type Riem 3 in the potential V are needed for having a renormalizable theory. 

To make the discussion of the redefinition theorem more transparent we present here an example of its 
use. 

Example. For the finite four-dimensional theory m, the following choice of F makes © tree-level 
equivalent to the Einstein-Hilbert theory (all the graviton scattering amplitudes are the same): 

e H(- n A ) _ i i e H(-a A ) _ i 

pm'.pa = g np g ™ -- -g^g pG ---+ s x g pv g pa RU{R-) + s 2 R p,y U{R p ' y ■) . (70) 

(When the action of the field redefinition (contraction F plJ ' pa R prj ) is evaluated, we substitute the center dot 
in the last two terms above by the Ricci tensor on the right R pa -) 

Remark. We here showed that S gr and Seh are tree-level equivalent and all the on-shell scattering 
amplitudes can be equivalently calculated using one or the other theory. However, off-shell amplitudes 
do not match, because in proving the theorem we made crucial use of the equations of motion Ric = 0, 
that uniquely characterize the perturbative graviton field in vacuum in both the theories. Similarly in the 
domain of classical field theory other non-perturbative solutions exist in S gT7 which are not shared by Seh- 
Therefore the two theories are equivalent only in the framework of perturbation theory, and full matching 
of the amplitudes happens only at tree-level. 

As reviewed in section H. 31 we can introduce other operators, local or nonlocal, to make the theory © 
finite at quantum level. In this section making use of the field redefinition theorem we proved that the 
killer operators in the action © do not give tree-level contributions to the n— graviton on-shell scattering 
amplitudes. However, © is not unique and we are free to introduce other killers that can affect the n-point 
functions. One example of such killers is the following quartic operator in the Riemann tensor, 

S 4 Riem 2 D 7 ~ 2 Riem 2 , (71) 

which gives contribution to the four and n— graviton scattering amplitudes (n > 4), as shown in the proof 
above, because we can extract from it at least four-graviton fields around flat spacetime (we remind that 
we get Riem ~ 0(h ), while Ric ~ 0(h 2 ) when expanding the metric around Minkowski spacetime with 
on-shell metric fluctuations). 

Beyond tree-level amplitudes. Finally, we expand about the operators we expect beyond the simplest 
tree-level computation. In this short subsection the main emphasis is placed on the finite terms and their 
contributions to the scattering amplitudes at loop levels. Of course it is known that such terms are not 
universal and the results for the amplitudes are unambiguous only if some renormalization conditions at 
some energy scale are fixed (below it is assumed that is is already done). 

At quantum level Einstein gravity is non-renormalizable and we expect contributions to n— graviton 
amplitudes from many other operators unlike the case of 1-loop super-renormalizable or finite theories. For 
the latter we expect to have an upper limit on the number of derivatives in the UV for the operators in 
quantum effective action (precisely 2q + 4 in D = 4), while for Einstein gravity we formally have up to 
infinite number of them. 

At one-loop a super-renormalizable theory in D = 4 gets extra nonlocal contributions that in the UV 
look like 

Rlog^-^-^R, Ric log Ric > (72) 
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and by virtue of the results in this paper they do not contribute to any graviton scattering amplitudes. 
Conversely the finite contribution, which arises in higher dimensions, 


Riem log 



Riem, 


(73) 


will give an explicit or implicit (after recasting in terms of the operators (1721) plus extra vertices) contribution. 
Moreover, we have an upper limit for the UV energy scaling of the scattering amplitudes evaluated with 
the quantum action. For the super-renormalizable theory without Riemann tensors, or a finite theory with 
killers not involving them, the upper limit is E 4 in D = 4, while in extra dimensions it is E D . 

With the field redefinition at hand we may attempt to compute the contribution even beyond tree-level 
in theories for which we can still apply our theorem. In this situation the two theories differ only by the 
Jacobian of the transformation that we need to include when we compute the partition function with the 
path integral method. As we proved above the action functionals are equivalent, i.e. S(g') = S'(g). In the 
path integral we also have to be careful about the change of the measure Jacobian of the transformation. 
The Jacobian of the field redefinition transformation g —> g’ is given by 


J = Det-^ = Det 


Sg^(x) 

s 9' p Ax') 


Tr log 


^9’pai.x') ' 


(74) 


Its divergent part vanishes in DIMREG scheme for the case of any analytic field redefinition. However, 
the finite non-analytic contributions to the quantum action, if any, take part in the Jacobian that we can 
compute perturbatively. For this purpose we introduce two ghost-like auxiliary fields (similar to Faddeev- 
Popov ghosts), which have fermionic statistics [34]. Next, we can consider Feynman loop diagrams with 
these fields and we compute the contribution of the Jacobian to the scattering amplitudes. At the zero-loop 
order (tree-level) the Jacobian is one and we should do not worry about the inclusion of diagrams with the 
new ghosts in the amplitudes. The contribution of the Jacobian starts at the one-loop order. 

In D = 5 we may expect the following finite terms in the effective action, 

qo R\/— □ R , q -2 RicV—□ Ric , q± Riem\/—□ Riem . (75) 


The non-analytic functions (like the above square root) of the covariant box operator appear due to dimen¬ 
sional reasons in any odd dimension. Only the last operator with two Riemann tensors will contribute to 
the amplitudes, and the coefficient q 4 will keep track of all the other form factors present in the classical or 
quantum action as explicitly evaluated in section |U 

We do not know at the moment the other finite contributions to the quantum action, but we know, as 
already mentioned, what the upper limit for the energy scaling of the finite contributions to the quantum 
action is: E 4 in D = 4. In the high energy regime (E 3> ~ Mp ~ A) we expect these finite contributions 

to approach the form of the following operators: 


1 1 9 
1Z—1Z—1Z 2 , 

□ □ 

111 , 

1Z—1Z—1Z—1Z 2 , 

□ □ □ 

... —iz 2 , 

□ □ □ □ 


(76) 


which probably contribute to the scattering amplitudes. 

Anyhow, only an explicit calculation of the one-loop scattering amplitudes will tell us the form of such 
corrections. 


6 . Conclusions 

In this paper we have performed a tree-level computation of on-shell four-graviton scattering amplitudes 
in the context of higher derivative Stelle theories and nonlocal gravitational theories quadratic in the curva¬ 
ture with non-locality specified by form factors. The theories in the first class are known to be generically 
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non-unitary due to the appearance of non-physical poles in the spectrum | 6 j. Conversely, the second class 
of theories, under some specific choice of the form factors, have been proven to be goo d candidates for a 
ghost-free and super-renormalizable or finite theory of quantum gravity 11. HEM 3E3 • I 11 both cases we 


have checked that, in the absence of terms quadratic in the Riemann tensor, the amplitudes coincide with 
the ones found in Einstein theory. Furthermore, the four-graviton scattering amplitudes in Weyl conformal 
gravity are identically zero. 

We provided an explanation of this result on the basis of a field redefinition map of quadratic gravity into 
the usual Einstein-Hilbert action [32| . This map is defined perturbatively at all orders in the parameters 
appearing in quadratic gravity. This ensures the two theories, sharing the same unperturbed action, are 
completely equivalent from the point of view of the tree-level amplitudes. More specifically this equivalence 
holds only as long as the two theories are supposed to have the same free spectrum. We cannot use the map 
to address the computation of observable quantities involving poles not present in Einstein theory. 

The idea that field redefinitions do not affect the physical S matrix was actually discussed in the context 
of quantum field theory, specifically renormalization theory, a long time ago and it actually lies at the core 
of the fundamental results about the quantum divergences of Einstein gravity 


the field redefinition very similar to the one in this paper can be found in 



An application of 
the authors argue 


that, just thanks to field redefinition, the order a'k 4 four-graviton amplitudes calculated in string theory 
can be derived from an effective action, where only the Gauss-Bonnet density shows up, consistently with 
the previous predictions in [56| . The field redefinition was also used to prove that in any higher derivative 
effective gravitational action the effective graviton propagator is always without ghosts. Similar results are 
derived at large in [33|. Our calculation is actually very near in the spirit to the one in [57J, where the 


effective quadratic action, reproducing the string theory amplitudes at the order a'k , is determined by only 
considering the 0(fc 4 ) amplitudes. However, in the approximation considered, the presence of the additional 
poles in the propagator contributes only linearly in 70 and 72 . 

In this paper wc have considered the full non-linear dependence including the full propagator of quadratic 
gravity and checked the redefinition theorem to all orders in 70 and 72 . When the Riemann square term con¬ 
tributes to the interaction vertices (in D > 4) we have found additional terms in the amplitudes, depending 
on 70 and 72 , which were neglected in previous computations. This is justified because the previous compu¬ 
tations were to order 0 (fc 4 ), while the first non-vanishing vertices’ contribution are possible to be derived 
from terms at least cubic in curvature (so at least of order 0(fc 6 )). Furthermore, the on-shell four-graviton 
amplitudes for a large class of weakly nonlocal gravity theories have been computed drawing on the fact 
that the form factors present in the action can be treated without much effort when on-shell gravitons are 
considered 
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We emphasize that our results differ from the standard ones obtained in Einstein gravity 
for the case of Gauss-Bonnet term in higher dimensions and also in the case of generalized Gauss-Bonnet 
terms (sandwiching a function of the d’Alembertian □ operator (1591) 1 in D > 4. In principle, this allows us 
to determine the form factors by comparing the results for amplitudes with hypothetical experimental data 
on graviton scattering. In particular we are able to determine the coefficients in the most general quadratic 
Stelle theory in D = 6 dimensions. Among other extensions of our results we mentioned amplitudes with 
any number of external gravitons (n > 4) and our expectations about results beyond the tree-level. 

Finally, although the original motivation of this study was to evaluate scattering amplitudes in a par¬ 
ticular class of weakly nonlocal theories of gravity, the outcome of the paper is a general feature of any 
higher derivative local or nonlocal gravitational theory: Einstein quantum gravity [331 ]. conformal gravity, 
effective string theory, and local or nonlocal higher derivative super-renormalizable theories j4f|. We are 
now motivated to introduce more advanced techniques in order to calculate n —point functions in local or 
nonlocal gravity involving the Riemann or Weyl tensors fbsl ] . 
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Appendix A. Details on the propagator 


In this appendix we will closely follow the procedure employed in [50(. Lorentz covariance and Bose symmetry 
allow us to rewrite the kinetic operator O a p n s(k) in (j5j) making use of the following basis, 


kp k(j 


Al — VpisVper 5 A 2 — ( f j jip'lisrr A '/ prr'/is p) j A 3 — 'I pis ,,, T Vper 


kpk v 


A 4 - 


k u k a kj/kp 

Vpp~Jj2 ^ Vpo 


k 2 


Vvp 


knkfj knkp 

— - L 


k 2 


k 2 


k 2 ’ 

kuki/k n k ( t 


k 2 


(k 2 Y 


(A. 1 ) 

(A.2) 


We note that O a/ 3 ^s(k) = /,A,; does not depend on 70 , 72 , and 74 , but only on the two linear 

combinations 70 — 74 and 72 + 474 , because the integral 


d°x {R^%-m^ 2 + R^ 2 ) 


(A.3) 


is identically zero in any dimension D. From action (I53D we read the coefficients 

/1 = (-2nn 2 k 2 ) i + 2fc 2 ^(70 - 74 ) + j (72 + 474 )^ , / 2 = (-2k^ 2 /c 2 ) ^-0 [l - fc 2 (72 + 4 74 )] , 


(70 - 74 ) + j (72 + 474 ) 


h = (-2 K~ D 2 k 2 ) (-2fc 2 ) 
h = (- 2 « 7 ) 2 fc 2 ) [ 2 (70 - 74 ) + (72 + 474 )] 


, / 4 = (- 2 K D 2 fc 2 ) (-/c 2 ) (72 + 474 ), 


(A.4) 


In order to define the graviton propagator we have to invert the kinetic operator O. These calculations are 


most conveniently carried out in terms of the Barnes-Rivers operators [5lL 
rank-two tensors. The complete set of D-dimensional operators is given by 


in the space of symmetric 


Ppis,pa — C) (Ypp^iscr + Op a U> V p + O^pLOpa + O va isJpp) j Ppv,pa — r, {9pp9isa + 9p<j9isp) \ 


piy^pa 7 


}O 0 


1 


-9,„,0 n 


} 0 c o 


jOOuj 


jOuO 


_ _ _ , , . . pu _ pyjuuj pu 

fils',per p _ 2 w fiv w per 1 J pi/,per ^per 5 1 p,u,per 1 pis',per ■ J pis',per 1 


-)§0oj 


1 


p 


,0^0 


1 


- 1 pis,per — ^ Vpissspcr , ^ P i/, P cr — ^ P<T ’ 

where and w pcr are the usual transverse and longitudinal vector projection operators 


0 /Jl/ = Vpis 


kpK 

“fc 2 ”’ 




kpk„ 

k 2 


that satisfy the relations: 0 pp 0 p = 6 p V , UppUj? = and 0 pp w p = 0. The operators P 1 , p 2 : P oe and P°“ 
are idempotent, mutually orthogonal, and satisfy the following completeness relation, 


\P + P + P +P ] fJtl/pcr — 2 ( T lpp r bscr + VpaVisp) — Ip 


They project out the spin-1, spin-2, and two spin-0 parts of the field. The two spin-0 transfer operators are 
such that 


pO pi _ plpO — p0p2 _ p2p0 _ q (pO) 2 _ p08 _|_ pOw 
pOu pO pO pOO pOtoO p06 pO pO pOcu pOOoj 

We decompose the operator 0 in the projectors basis 

O a p„s(Je) = ci P 1 + c 2 P 2 + c£P°“ + c®P oe + c 0 P°, 


(A.5) 

(A. 6 ) 

(A.7) 
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where the coefficients ci, C 2 , Cg , eg, Co are obtained using the tensorial identities 

x 1 = (D- 1) p 08 + p 0ui + Vd~Tp° , X 2 = P 1 + P 2 + P 08 + P 0ui 
x 3 = Vd~1P° + 2 P 0u , X 4 = ip 1 + P°“ , A 5 = P°“ . (A. 8 ) 

The coefficients <7 are explicitly: 

ci = h + = (- 2 Ko 2 k 2 ) i c 2 = f 2 = 2 K" 2 fc 2 (1 - k 2 (72 + 474 )) , 


eg = (0 - 1) /1 + h = {~^D 2 k 2 ) - [(D - 3) + 2 k 2 (4 (D - 1) (70 - 74 ) + D (72 + 4 74 ))] , 


(A.9) 


Cq — /l + fl + 2/3 + /4 + /5 — ( — 2 K D 2 k 2 ) ( — — 


Co = VD ^ l ( f 1 + f 3 ) = (- 2 K ~ D 2 k 2 ) - s / D ~= 1 . 


The problem of finding O a p s (k) = s 4 P l + s 2 P 2 + SqP 0ui + s 8 P 08 + sqP° boils down to solving the linear 
system 


o-o - 1 = 
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Using the echelon matrix form [5 
/ci 0 0 0 0 1 \ 
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the propagator is given by 
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Appendix B. Useful expansions 

We list the expansions that are necessary to determine the propagator and vertices used in the paper. The 
expansion of the metric around flat spacetime is defined as 

H/ii/ — A (U.l) 

We will always assume as an on-shell field satisfying the conditions Oh^ = 0, = 0, = 0. 

This is convenient to reduce the expressions to a compact form which is the one we actually need in most 
computations for on-shell gravitons amplitudes. When we want to refer to unconstrained off-shell fields we 
will adopt the notation h . 

g ^(i) = _ h iu , ; 

/"( 2) = h> M h a V , 

(V^) (2) = ~lh ab h ab , 

P (1) = d a d b h ab - UI& , 

P (2) = -d b h ac d c h ab + \d c h ab d c h ab , 
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i? (4) = -12h ab h cd d b h de d c h a e + 18 h a c h ab d b h de d c h de + 24 h a c h ab h de d c d b h de - 24 h a c h ab h de d c d e h bd 
+36 h ab h cd d c h a e d d h be - 24 h ab h cd d b h a e d d h ce - 24 h a c h ab h de d e d c h bd + 24 h a c h ab h de d e d d h bc 
+24 h ab h cd d d h ce 8 e h ab - Qh ab h cd d e h cd d e h ab - 24 h ab h cd 8 d h be d e h ac + 18 h ab h cd d e h bd d e h ac 
-24 h a c h ab d c h de d e h b d - 12 h a c h ab d d h ce d e h b d + 3 6h a c h ab d e h cd d e h b d , 

R$ = \ + d p d a h av + dvPhvr - Oh, J , 

R$ = \d p h ab d u h ab + h ab ( d b d a h pv + d p d v h ab - d b d p h ua ~ d b d v h pa ) + d b h p a (, d b h va - d a h ub ) , 


— 2 ( phvi 7 + dfi&ahup dvdphj lrj d u d a h pp ) , 


(B. 2 ) 


R$ga = -\d a K a d a h^p + \d a K p d a hp a + \d a h vrj d p h pa - \d a Kpdph aa 
-~\d a hp<jd v hpa + \dph aa d v h p a + \d a hp P d u h aa - \dp,h p a d u h cja 
h f) h a — —f)h f) h a — -f) a h f) h -f) h r) h a 

i 2 w a'2 L/l ' t<TaU P ''L 1 2 U ,6 M crL/ P /t ' iya ' 2 L/ I- L 0(7(1 

-\d a K p d a hp a + \d v h pa d G hp a + \d p h va d G hp a 

+\d a hppd a K a - \d p hp a d a h ua - \dph pa d a K a , (B.3) 

= -§ h ab dahp<jd b Kp + | h ab d a h IJ . p d b h„ tT - | h ab d a K a dph pb + | h ab d a h up dph ah 
h ab dphcr b d u hpa + | h ab d a hp a d v h pb - \h ab d a hp P d v h ab + | h ab dp,h pa d v h ab 
-§ h ab d b h va d p hp, a + \h ab dvh Gb dphp la - \h ab d p h ab d p h ua + | h ab d a h pa d p h vb 
+ \h ab d b h V pd a h pa - \h ab d v hp b d a h pa - | h ab d p K b d a h pa 

+| h ab dph pb d a K a - f h ab d a hp P d a h„ b + | h ab d p h pa d G h vb . (B.4) 
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